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On the Numerical Solutiofi of Integral- Equations, 

Bj E. T, Whittaker, F.R.S. 

(Received December 27, 1917.) 

1 1. Introductory, 
The present communication is concerned with integral-equations of Abel's 

I (l>(s)fc{x—s)ds =f{x), (1) 

Jo 

and of Poisson's type 

where k {x) is a given function called the nudetis^ f{ce) is also a given function, 
and <f>(x) is the unknown function which is to be determined. The object of 
the work is to obtain solutions of these equations in forms which can be 
made the basis of numerical calculation. 

Theoretical solutions of both these equations, in the form of infinite series, 
are well known, and have been fully discussed by Volterra* and others. 
But in these solutions the 7ith term of the series is a multiple integral 
involving (n — 1) integrations with variable limits : and although such series 
are valuable for the light they throw on the general properties of the 
solution, it is obvious that they cannot, except in very special cases, be used 
in order to compute values of the solution numerically. The only case, so 
far as I am aware, in which a solution of an equation of one of the above 
types has been obtained in a form adapted for practical ends is AbeFs 
original special form of equation (1), 

^ <!> (s) ds 



■■/(x), HXiXl, /(0) = 0), 
for which he gavef the solution 



o(r^-s)P 



1 . 



X 



f (5) ^^ 



■-V 



•^- ^Torino Atti,' vol. 31, pp. 311, 400, 557, 693 (1896). For fuller references cf. 
H. Bateman, " Report on the Theory of Integral Equations," 'Brit. Assoc. Report,' 1910. 

t ^ CEuvres,' (ed. 1881), p. 11 (1823) and p. 97 (1826). The fundamental meaning of 
AbeFs result is most clearly seen if the integrals whicli occur in it are interpreted as in 
the theory of generalised differentiation : if i|r(^;) is written for V{\.~'p)^{pG)^ Abel's 
formula reduces to the simple statement that if 

then ^{x)=^{~-\'^f'{x), 
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Whenf(x) is given, the values oif\s) and of the integral last written may 
be obtained without difficulty by the ordinary processes of interpolation and 
numerical integration.* 

Kecently, two friends, one a seismologist and the other an actuary, have 
enquired of me whether the integral-equations (1) and (2), which had 
occurred in their reseeirches, could be solved in such a way as to obtain 
numerical results when the functions /c (x) and f(o^) are known (tabulated) 
functions. It was under the stimulus of these enquiries that the methods of 
solution which occupy the following pages were devised. It will be seen 
that I have departed altogether from the customary methods of solution by 
infinite series whose terms are multiple integrals, and on this account the 
new solutions, which are formulated in Theorems 1--5 below, and by which 
the unknown function may be determined numerically, may perhaps be 
found to be not without interest from the standpoint of pure theory. 

§ 2. Solution of Integral- Equations of AbeVs Type. 
Considering first the generalised Abel's integral-equation, 



I 



X 

(f> (S) fc{x~-s) ds := / (X), (1) 





we need only consider the case when the nucleus k (x) becomes infinite at 
X = 0, for in the simpler case when the nucleus is finite at x = 0, the 
equation (1) may be reduced immediately, as we shall see later, to the 
type (2), and so may be dealt with by the methods which are given subse- 
quently in the paper. 

We shall, then, suppose fc (x) to be such that x^/c (x) is finite and not zero 
at X = 0, where ;p lies between and 1. 

Now if the nucleus /c (x), which is supposed to be given by a numerical 
Table, has this character, so that it becomes infinite like x"^ at x = 0, but is 
finite for other values of x within the range of integration, we can in general 
(by use of ISTewton's or some other interpolation-formula) represent the 
function xP/c{x) over the range in question by a polynomial in x; the degree 
of this polynomial will depend on the nature of fc (x) and the order of 
accuracy to which the work is to be carried. We may, then, assume for the 
nucleus k (x) an analytical expression of the form 

fc (x) = x"^ (ao-t chx -f a2X^ -f . . . -f a^x""), (0 < p < 1). (3) 

'■^' It is possible that the elegant solution in power-series which was given by Sonine, 
' Acta Math.j' vol. 4, p, 171 (1884), for his generalised form of Abel's integral-equation,, 
may also be ntilised for numerical calculation. 
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Now assume* a solution of the form 



^ (x) = f {s) ^{x—s) dsy 



(4) 







where K {x) is called the solving function. We have on eliminating (^ (x 
between this equation and (1) 

f(x) = f/c (x-s) I fV (0 K (s-t) dt~\ds, 
Jo L Jo J 

or, inverting the order of integration, 

f (t) dt =: \f(t)\ K {X-^S) K {S^t) ds ydt, 

Jo L Jt J 

Since /^ (i^) is an arbitrary function, this gives 



i? 



fc(x—s)lL(S'—t)ds= 1. 



Writing s = t-^it, this becomes 



*x—t 







or 



/c{x---t — u)K (u) du ■■= 1, 

•a 

K (a—u) K (?^) (^'Zi = 1, 





where a is arbitrary. From this equation the solving function K (x) is to be 
determined. 

By (3) this may be written 

(a—uyp {ao-h ai(a—u) + chia -nf -}-,,. -i- ttnia—tiy} K(ti)du = 1. (5) 

Jo ' 



JO 



IsFow let 



'U 







K(u)du be denoted by Ki(u); let Ki (^^) ^-^z. be denoted by 



'u 



K2(u), and so on. Then integrating by parts we have 

'a /•« 

(a—ity--PK(tc)du = (l—p) (a—uypKiiu)^^, 

o' Jo 



Jo 



*a 



(a~uy-PK(Mldul=^ {2--'p) {a~uf--PK^{u)du, 
Jo 



'a 



Jo 
and so on. Thus equation (5) may be written 

(^a - ^l)-P { ao K (to) -f (1 -^) ^iKi (^) + ( 1 -p) (2 -~jp) asKs (t^) -f . . . 
Jo 

+ (1 —-p) (2 — jt>). . . (?i --p) an Kn (u) } du = 1 . 
"^ The legitimacy of this may be inferred at once from Volterra's theory. 

2 G 2 



(6) 
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*a 



Now if I {a--u) ^h{u)du = 1, where a is arbitrary and h{u) does not 
Jo 

involve a, we have by writing % = as 



aP-'^ = I (l — s)-Ph(as)ds, 
Jo , 



and, since a is arbitrary, this shows that h (u) = Ct^?'""^ where C is independent 
oiu. 

Substituting in the last equation, we have 



\\i-sy 

Jo 



sm pir 
and therefore 

TT 

Applying this result to equation (6), we have 

ao K (i^) + (1 —p) ai Ki (ic) -\- (1 —p)(2 —p) a^ K2 {n) + . . . 

IT 

or, writing y (ti) for K„ (ti) 

+ {l-p)(2-p)...(n-p) a„y ^ -'-Hi£!I . uP'-K (7) 

This is a linear differential equation in y, and K^ (?0 is that solution of it 
which vanishes, together with its first (n—l) differential coefficients, when 
■I'u vanishes. 

Now let the polynomial 

aQX'' + (1 — p) ai9f~^-{-(i —p) (2 —p) a2X"~^ + . . . + (1 — j/) (2 —p), .,{n —p) an 

be denoted by 'F{x): and let its n roots (supposed for the present to be 
distinct) be oc, /B, 7, . . . , p» Tiien it is known from the general theory of 
linear differential equations that the solution of (7) which vanishes, together 
with its first (n — 1) differential coefficients, when ti vanishes is 

^aX n^X gvX 

where h (,«) = ^^^ + ^^^^ + . . • + j^^. (9) 
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From (8) we have by successive differentiation 
Kn-i (u) = ^^5^ fV-1 A' (ti--t) dt, since h (0) = 0, 







Kn-2(u) = 5i^^^ tP-^^h'' (u---t)dt, since h' (G) = 0, 

TT Jo 



Jo 



TT 



(^,) := 5Hi^/!?^_!+ (ViA<^0 0^-0 ^^T. since A<«-i)(0) = -/** 

77" L ^0 Jo J ^0 

Substituting for h{x) from (9) we have K (i^) = (I/tt) sin^7rL{a?), where 

/j,jp— 1 .,?^ Cx On rx 

-L{x)=.'^ + i5/7-T6"M tP-^e--'dt + r^r-^e^^\ p-^e-^dt-^r... 

% JB (a) Jo ' F O) Jo 

1 W Jo 
rx rax 

.Now tP-^e-'^UU = a-P s^-^e-'ds; 

Jo Jo 

and thus if the function e^' sP~'^e~Uls 



1?"^' 



(which ia well known under the name of the Incomjplete Gam/ma- FuticUoki) 
be denoted by 7^ {x), we have 

Combining our results, we have 

Theorem 1. — The sohction of the integral-equation 

fX 

^ (s) k(X'— s) ds = / (x) 

where the nucletts fc{x) is supposed to he given minfierically and to have been 
expressed ly the ordinary methods of interpolation in the form 

n (pc) = dcT'P (ao + a\X + a^o^ + . . . + <x?^^*0> (^ < i^ <C 1) 
is 6 (x) = ^^^-^ f/ (s) L (^-s) ds, (10) 

TT Jo 

'where 

^ For A(?i-i) (0) = simi of residues of t'^-'^l¥{t) at its singularities «, 0, ..., z/ 

^ - (residue of this function at oo ) == I/^q. 
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unci ivhere a, ^^ ,.,, v, are the roots of the algehraic equation 
F {x) = oo^^ -f (1 — p) aix^--^^ (1 — p) (2 ~-^) a2;r^"2 + . . . 

+ (l-.jp)(2--^)...(^-~^)6^^ = 0, (12) 
ant^ 7j3 (x) denotes the incomplete Gamma-fimction 

sP'~^e~^ ds, 



This may be regarded as a direct extension of Abel's original formula, 
which may be derived from it by taking n ~ 0. It expresses the solution 
of the integral-equation in a finite form in terms of the incomplete 
r~function. ■ The incomplete F-f unctions which occur in the solution all 
have the same parameter jp, and are, therefore^ really all the same function, 
with different values of the argument ; its values may be tabulated from any 
of the expansions which were given for it by Legendre,* such as the 
absolutely convergent expansion 

p jp(p + l) jp(p + l)(p + 2) 
or (for large positive values of oS) the asymptotic expansion 

Whqn this function has been tabulated, and the algebraic equation (12) 
solved,f the function L(^) may readily be tabulated fl'om equation (11), and 
then the required function ^{x) is given at once by quadrature from 
equation (10). In this way Theorem 1 yields a numerical solution of the 
generalised Abel's equation. 

It is obvious that when the polynomial F(^) has a pair of equal roots, 
certain modifications must be made in the above solution, but it does not 
seem necessary to set these forth in detail here. 

It may be remarked that by making the degree n of the polynomial in (3) 
increase indefinitely, we obtain in the limit the solving-function of the 
integral equation, with an arbitrary function as nucleus, in the form of an 
infinite series of incomplete F-functions of the same parameter. In con- 
nection with such a solution it would be necessary to discuss convergence, 
^tc, and it is not proposed to undertake this in the present paper ; but the 
matter seems worthy to be mentioned in passing, as series of incomplete 
F-functions have not (so far as I know) presented themselves hitherto in 

■^- 'Exerc. de Calc. Int./ vol. 1, pp. 338-343 (1811). 

t For this Newton's method will probably be found in most cases the most convenient, 
if n>Z ', or if some of the roots are complex, the Lobachevsky-Graeffe method. 
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analysis. They may evidently be regarded as an extension of Dirichlet's 
series, i.e.^ series of the type 

which have received much attention in recent years. 

§ 3. Solution of Equations of Poisson's Type. 

When the nucleus k {x) of equation (1), instead of being infinite at ^ = as 
we have supposed hitherto, possesses a finite differential coefficient over the 
range of integration, we have by differentiating equation (1) 

^ {x) /c (0) + I (^ (s\k {x—s) ds = f {x) 

which is an integral equation included in Poisson's type (2).* We shall, 
therefore, now pass on to consider integral-equations of Poisson's type, which 
we shall take in the form 

^{x)-\- <p (s) fc(x--s) ds = f{x). (2) 

Jo . 

Since the nucleus /c{x) is supposed to be specified by a table of finite 
numerical values over the range of values of x considered, we may apply 
Prony's method of interpolation by exponentials in order to represent it 
analytically in the form of a sum of /x exponentials 

/c (x) =• FeP"^ + Q^^'' + Re''^ + . . . + V^''^ (13) 

where (P, Q, R, ..., V, ^, q, r, ...,v) are constants which are chosen so as 
to give the closest possible representation of the given numerical values. 

Although Prony's method is more than a century old, it does not appear to be widely 
known or to have found its way into any text-book ; and, as his original paper is perhaps 
not accessible to many readers, I may be justified in giving here a brief notice of it. 

Suppose that k {cc) is given numerically for a certain range of values of x. Take any 
set of values of x equally spaced within this range, say ^ = 0, ca, Sw, 3a), 4ci), ..., and let 
the corresponding values of k{x) be k^, kj, ^2,^3, .... Now if k{x) could be represented 
exactly in the form of a sum of fx exponentials, say, 

Pe*^^ + Qe«^ -f B.e''' + . . . + Ve^^, 

then K (x) would satisfy a linear difference-equation of the form 

AKu+/x + Bi<ri+^~i-{-C.«C7i+^-24-. 4-Mk^ = 0, 

where the roots of the algebraic equation 

Axf^-\-'Bxi^~'^-hCxy-~''^-{' + M = 

would be 

(?p% ee^^ grw^ qv<o^ 



* It need scarcely be said that if k(0) vanishes we differentiate again. 
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Prony's method, which is based on this fact, is to write down a set of linear equations^ 

Ak^ +B#c^-i + C^,u-2+ + M^0 = 

Ak/j + i + Bku +C#£/a-i4-......+M/Ci = 

A%+2 4-Bk^+i + Ck^ +....., +Mk2 = 



where the quantities Kq, Kj, ^2, % ..., are known, since k{cg) is a known tabulated 
function, and by the ordinary method of Least Squares to find the values of 
A, B, C, ..., M, which most nearly satisfy the equations; then with these values of 
A, B, C, ..., M, to form the algebraic equation 

Acey- +B.i7M-i-|-C^M--2+...4-M =- 0, 

and find its roots ; these roots will be ev^^ eQ^, e''^^ ... <3^"% and thus p, q, r, ..., 'y, are 
determined. Knowing p, q, /•, ..., v, we have a set of linear equations to determine 
the coefficients P, Q, K, ..., V, and these also are to be solved by the method of Least 
Squares. 

Taking then this form (13) for the nucleus k (x), we shall show that the 
integral-equation (2) may be satisfied by a solution of the form 

(j)(oi) = f{x)— K(fc— s) f (s) (Is, (14) 

Jo 

where the solving function K (x) is also a sum of /a exponentials, say 

K {x) = A^^'^ + B^^-^ + C^Y^ 4- . . . + N'6''^'. (15) 

To prove this, we remark first that the existence- theorems established by 
Volterra justify us in assuming for the solution the form (14), where K {x) is 
now the function to be determined. 

In (14) put fc {x) ior f(x) : thus 

(j) (jjo) = f{;{x)— \ 'K(x—'S)h: (s) ds, 

Jo 

which gives the value of (^(c«) corresponding to this value oif(x). 
Putting (x—s) for s in the integral, we have 

(f}(fc) z=z fc-{x)-~ IL {s) K {x — s) ds. 

Jo 

Comparing this with the integral-equation (2), after replacing f{x?) by Kix) 

in the latter, we have 

^ (cc) = K {x)^ 

and therefore the pair of functions 

^ (,7j) = K (cc), f{x) - K (x\ 

satisfy the integral-equation : that is to say, 

.K(i^)+ lL{s)K{X'~s)ds = fc{x). (16) 
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In this equation substitute the vahie (13) for fc{x) and the value (15) for 
K (m). Thus we have 

Jo 
Equating coefiieients of e"^' on the two sides of this equation, we have 



^ +_!!_ + ... + -JL_+ 1 = 0. 



«_p «-? 



a — v 



Similarly by equating coefficients of e^^ 



p n \T 

-|- __ — __». — |— ^ ^ ^ — 1« „_«^ — , — . ~^ j_ — ^j^ 



^-P ^-f 



^ — '-y 



and so on: these equations show that a, 0, % ..,, p, are the roots of the 
algebraic equation in x 






V 



»;— p 33— g X — r 



iAj ""*" 'V 






(17) 



This enables us to determine a, 0, y, ,,,, 

Next equating coefficients of e^^' on the two sides of the equation, we have 



and similarly 



I -»-^ I ^-^ t I XN , -| /\ 



a 



-P ^-P 7-P 



y-P 






*-»*«•»» ^« '»»****9**m**0»**9»ikmm***«#m%»tfm *»««•»# 



AB.C, ,N,T 



v. 



a 



•V ^-— -?; fy-~^ 



V — V 



J 



(18) 



Since («, ^6, 7, ..., 2/) and (j?, r^, r, ..., i?) are known, these equations (18) 
enable us to determine A, B, 0, ..., ST, and we see that if the constants 
(a, ^, 7, ..., 2;) and (A, B, C, ..., N) are determined by equations (17) 
and (18), the equation (16) is satisfied by the value (15) of "K{m\ 

The value of K(a?) may be obtained in a more explicit form in the 
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following manner. If we eliminate A, B, 0, ..., N, determinantly from 
the equations (15) and (18), we have 



1 



r 



e 



a -—p /3 —p V -—p 

111 



oc-q^-q v-q 



• «• ••4»4' *»*9%^*'*9* SO 



^••VilVtOi* 



1 






a 



•V 0—v ' r— ?; 



1 



1 



1 1 

0-q"p-q 



» a a » » i 



9«4"»9»'«»tt«a«a«« 



1 



1 



'€ 



^x 



-I- • * « ' 



^-F 



V--P 



1^1 



^-9. ^-9. 



• ••« • ••• ••«• •••» 



• «••«•< 



• • • • 9 i 



_*L « » 9 ' ' ' 



a — v 



y^l) 



^■—V Iz — V 

The determinants which occur in this equation are of the kind known as 
alternants, and may be faotorised by known methods.^* Performing the 
factorisation, we have 



Jul ( «^ J 



(ct-p(a-^q)(a-r)...(oc'~'V) ^^^ 



e 



((3d — 18) (ct — 7 j. . .((5t — f) 



,vx 



_ (2/ -— p) (p--q)(v—r)..,{p'- v) 

Combining our results we have 

Theorem 2, — The solution of the integral-equation 

<f> (5/) + j <^ (s) k: (x —s)ds = / {x), 
Jo 

lohere the 7iudeus ic(m) is supposed to he given numerically, and to ham been 
expressed hy Pron/ys method of interpolation in the form 

Cx 

(f)(x) =^ /(x)— I K(x—s)f(s)ds, 



where 



K(m) 



and where a, ^, 7,.-, v, <^^'^" the foots of the algebraic equation in m 

P 



e 



vx 



(19) 



- + 



E V 



x—p X — q X — r X — v 

* The evaluation of alternants of this type is due to Cauchy, 'Exercices d'An.,' vol. 2, 
p. 151 (1841). 
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It is obvious from the last equation that if p, q^r^ ..., v, are arranged in 
ascending order of magnitude, and, if P, Q, E, ..., V, are all positive, then the 
lowest root a is less than p, the next root ^ is between p and g, and so on. 

If the number of exponential terms in (13) is supposed to increase 
indefinitely, the representation of the nucleus k{x) becomes a Diriohlet's 
series, and by (19) the solving function K {x) is then also a Diriohlet's 
series, formed with exponentials e"^, e^^, ^y^, ..., whose exponents oc,/3, 7, ...,., 
are the roots of equation (17), which now becomes a transcendental equation, 
A rigorous examination of convergence, etc., in this limit-process would be 
necessary to establish the theorem which appears to be indicated, namely, that 
in the solution of a Poisson's integral-equation whose nucleus is expressible as a 
Birichlefs series, the solving -function is also expressible as a Birichlefs series, but 
with a different set of exponents for the exponentials, 

§ 4. An Alternative Sohition of Integral-JEquations of Poisson's Type. 

Theorem 2 above supplies what I think will be found to be in general the 
most convenient method of solving integral-equations of Poisson's type 
numerically. But in certain cases the nucleus tcix) may be given as a 
polynomial, or it may happen, when /c {x) is given in the form of a numerical 
table, that it is preferred for some reason to apply ordinary interpolation 
and express ic{x) approximately as a polynomial, rather than to apply 
Prony's method of interpolation and express k {x) as a sum of exponentials.* 
We shall therefore now consider the problem of integrating an integral- 
equation of Poisson's type 

(f> {x) + ^{s)ic {x—s) ds =f (x) (2) 

Jo 

when the nucleus k(x) is expressed as a polynomial. 

The solution of this problem may be deduced as a limiting case from that 

solved in | 3. For, in equation (13), suppose that p, q, r, ..,, v (being n in 

number), each tends to zero, while P,Q, E, ..., V, increase indefinitely in such 

a way that 

P + Q + E+...+V == fco 

Vp + Q.$ + Er + . . . + Yv = /ci 



P^^-i + Qq"""^ + Er^""^ -h . . . + Yv""-^ = /c 



>- (20 



w- 



^ 



* One great advantage of Prony's method is that each exponential term involves two 
disposable constant^ (e.g., the term Pe^^ involves P and p\ whereas a term of a poly- 
nomial involves only one disposable constant (e.g., px^ involves only p), and therefore 
it is in general possible to obtain as high a degree of accuracy in approximating with n 
exponential terms as with a polynomial of 2^ terms. 
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where, /fo, ki, k^, ..., K„-i, are finite. Then from equation (13) we have 



tc (x) == PM -f'^^j-f c^-f. ... ) -f Q (l^qx-\- 



+ Y ( 1 + '^^^ + 






n^ .... I T^ . . . 



\ 2! • 

the terms in higher powers of x vanishing. 

Moreover, the equation (17), whose roots are a, ^, j, ..., v, may be written 

\tAj dj Jj / ■ tA/ i.{j tX/ / \tAy i/^ oO / 



which by (20) becomes 



/<^o , ^1 , ^^2 , , f<^n-l 



i 5 *T -o "T" • • • "T" T -»- • — ^; 



fAy lAy w t'Cy 



the terms in higher powers of {Ifx) vanishing, so that a, 0, j, ..., v, are now 
the roots of the algebraic equation 

The equation (19) for the determination of the solving function K {x) now 
becomes 

^^^ {oL-^^){a~-^),,.{of~v)^ 0-(x)0-7)...0-~F) 






and thus, collecting our results, we have 

Theorem 3. — The sohdion of the integral-equaUo^i 



"X 



^ (x) + ^ (s) f€{x— s) ds = / (x), 

tvhere the nucleus fc (x) is siqoposed to le given wiimieriecdlg, ami to ham been 
expressed hy the ordinary methods of interpolation in the form 

is j>{x)^f{x)-^\lL{x^s)f{s)ds, (22) 

where 

-^ W - (a^^)(a~~ry)...(^^Z;) (/i^-^ a) (^-7). . .O-^f) 



^n 



(f — a) (f — ^). . .{v — iJb) 



e-^, (23) 



On the Numerical Solution of IntegraUEquations . 379 

•mid whefG ot, j3, 7,..., v, are the roots of the algehra/ic equation in x 

£>/* + a:o^^"^4-a:i^"~^+... + a:«~i = 0. (24) 

If this equation (24) has a pair of equal roots, terms of the type xe^^ will 
occur in K {x\ thus, if 

K {x) = — 2p -{-p^x, 

we find K {x) = — j> (2 -{-px) e^"". 

An interesting expansion is obtained by making n increase indefinitely in 
this theorem. We then have the solving problem K(a?) of the general 
Poisson's integral-equation expressed in the form of a Dirichlet's series, or, 
at any rate, a series of exponentials with real or complex arguments. This 
appears to indicate a new field of analysis, in which Dirichlet's series present 
themselves naturally. But a thorough investigation of convergence would 
be necessary for the rigorous establishment of this result, and, indeed, we 
can show by simple examples that the expression obtained by making n tend 
to infinity will not necessarily be a Dirichlet^s series, even though, for all 
finite values of 71, it is a sum of exponentials. Eor instance, if 



Ti T r) Of r) 



m that ^0 = P> '^i = P^> a:2 = p% . . ., /Cn-i = ^:>^ 
we see from (23) and (24) that K (x) is equal to the sum of the residues of 

__ (l~p/Of>' 

1 — ^«+7^^+i 

at its n poles, which are the (n-i- l)th roots of unity other than unity itself. 
This sum of residues is evidently a sum of exponentials in x, one corre- 
sponding to each pole, so long as n is finite, but, when n increases indefinitely, 
we have 

.K{x) =z coefficient of 1/t in — (1— j?/^)e*^ = p, 

so that when the nucleus is k(x) = pe^"", the solving function is K(x) ~p, 
which is not a series of exponentials. 

§ 6. ^ Further AUernatim Solution of Integral-Equations of Foisson's Type, 

Hitherto we have supposed the nucleus of the integral-equation to be given 
by a numerical table, and to be represented, by use of the methods of the 
interpolation theory, as a sum of exponentials, or as a polynomial. If, 
however, in an equation of Poisson's type 

^(0?)+ fc{x—s)^(s)d$=f(x) 
Jo 
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the expansion of the nucleus k {x) as a Taylor's series is known^ say 



(25) 



then we shall show that the solving-function K(a?) may be written down at 
once as a Taylor's series. For, the solution of the integral-equation in terms 
of the solTing-function being 



rx 



^(x)=f{x)-' K (x— s)f (s) ds, 







we have already in equation (16) shown that 



rx 



K(x)+ K (s) K (x—s) ds = fc (x). 

Jo 

Putting X = in this, we have 

Ko == fcoy (26) 

where Ko, Ki, K2, ..., are used in order to denote the values of K{x) and its 
successive differential coefficients at :^ = 0. 
Differentiating (16), we have 

K' (it^-f KO K (x)+ Tk (s) k' (x-s) ds =: k' (x), (27) 

Putting ^ = in (27) 

Ki + /<:oKo == ^1. (28) 

Differentiating (27) 

K''{x)~\-fco'K'(x)-{-fciK(x)+ {\(s)K''(x--s)ds= ^'' (x). (29) 



Putting 02 = in (29), 

K2 + /^oKi + a:iKo = /f2. (30) 

If now the linear equations (26), (28), (30), ..., be solved for Ko^Ki^ K2, ..., 
we obtain 



^0 - ' 


Ko 








s:i- 






fCo 


1 








Kl 


Ko 


l2 = 


Kq 


1 







^1 


fCo 


1 






^2 


m 


Ko 



etc. 



Let us now consider the convergence of the series 



This series will evidently be absolutely convergent for all finite values of x,. 

provided the series 
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converges within a circle of non-zero radius. But the equations (26), (28), 
(30),..., may be comprehended in the single formal equation 

(/co + /ci^ + a:2^'^+ ...)(l + /^o« + /^i^^ + /^2«^+ ...) ~ Ko4-Ki^-|-K2i^2+ ... 
or 



Ko + Ki^^ + K2'«^ + Ks^s ^ ^ ^ ^ 



1 



X 



v_ 



1 + K(fiO + K\0? 4- IC<^^ + . . . 



and the expression on the right-hand side of this equation represents a 
holomorpliic function within a circle of non-zero radius having the origin as 
centre, provided the power-series 

converges within a circle of non-zero radius and so represents a function 
whose singularities and zeroes are all at a finite distance from the origin. 
Subject to this last condition, then (which is, as a matter of fact, unnecessarily 
stringent), the series 



X' 



X? 



converges absolutely for all finite values of x. We shall assume that it 
converges, since the computer will not make use of this method unless the 
convergence of this series is so rapid as to be obvious. Then, combining our 
results, we have 

Theorem 4. — The solution of the integral-equation 

'X 

^(x)-\- ^(s)K{x—s)ds=f(x) 
Jo 

where the nucleus k (x) is supposed to he expansible in a Taylor series 

k(x) = Ko-i-KiX-h^ +'^'^... 



%s 
where 



cj,(x) =fix)-^i\(x-^s)f(s)ds 



K(x) = /Co 



/Co 1 



X + 



Ko 1 
Ki Kq 1 
K2 K\ Ko 



0^ 

21 



Ko 1 

Ki Ko 1 

K2 Ki Kq 1 

Kz K2 Ki Ko 



,3 



X 

3l 



+ ... (31) 



For the benefit of those who may find it convenient to use this solution, it may be 
well to add some remarks on the computation of the determinantal coefficients which 
occur in it, A numerical determiuant should never, or scarcely ever, be evaluated by 
expanding it; it should be evaluated by reducing it successively to determinants of 
lower order, without expanding. To do this, we first notice whether any of the 
elements in the determinant is unity ; if not, we reduce one of the elements to unity 
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by dividing the row or colunin wMcli contains it by that element, 
to evahiate^ e.g.^ the fifth -order determinant 



Having done this, 



1 


(-h 


C?3 


a^ 


% 


h 


h 
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^5 


^1 


C2 


% 


H 


% 


di 


d^ 


^ 


d\ 


d. 


H 


09 


% 


% 


% 



the rule is : Strike out the row and column which contain the unit element, and 
subtract from each of the other elements the product of the elements which are 
situated at the feet of the perpendiculars fi^om this element on the deleted row and 
column. Thus the above determinant becomes 



h^—a^i h^-GL'^i b^--aj)i h^-UrJbi 

%-«2Cl %-«3^1 ^4-^4% %-%Cl 

e^~a^^ ^a-«3^i 64-^461 Or^-a^e^ 



and in the same way we may reduce this fourth-order determinant to a third-ordei^ 
determinant, and so on. 

If the unit element is oot the leading element, we must multiply the whole deter- 
minant by +1, according as the sum of the row-number of the deleted row and the 
column-number of the deleted column is even or odd. 



I 6. A Combmation of the Solution qf | 4 luith that of | 5. 

A form of solution which in some cases proves useful is obtained by com- 
bining the solution of § 4 with that of § 5. This happens when the graph of 
the nucleus k{o(:) in the part of it over which the integration takes place, is 
not very different from that of a polynomial of low degree, so that the first 
few^ of the coefficients /co, ^1, /^2, -..; are of preponderant importance as com- 
pared with those that succeed them. In this case it is advantageous to take 
out of the series (31) the terms which depend solely on these important 
coefficients, and, by summing them, to obtain a new form for 'K(x) which can 
be more readily computed ; this, as we shall now show, may be done, the 
extracted part of K(a)) being in the exponential form which was obtained 
in I 4. 

Suppose, for instance, that /cq and ici are important, but the succeeding 
coefficients /^2; a^s,-. are comparatively smalL The terms in (31) which 
depend only on kq and ici are evidently 



K\) 



Ko 1 






Ko 


1 





9 


Kl 
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"1 


21 
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i€i^ \ 
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K.i 


^0 
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fCi 


/Co 
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Ki 


Ki) 



<> 

X"-' 



3! 



'WB^IHHK 
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Kow if Ur, denote the determinant 



(-) 



n~\ 



Ko 


1 








0... 


fC\ 


%. 


1 





0... 





Ki 


fCii 


1 


0... 








K^l 


A^O 


.JL « • • 





{n rows) 
we have by expanding in terms of the elements of the first row 

The solution of this difference-equation is 

Un = A«^ -f- B/3^* 
where a and /3 are the roots of the quadratic ocP -{- fcox -h f€i =.0, and Aand B 
are independent of n. Moreover, since 



lOi 



fCQ = — (a-f-/S) = — 



and 

we see that 

and therefore 






a 



3__^8 



«-^' 



A = 



u^^ — 



Ai + 1 Qn-\-l 



a' 



The expression (32) may thus be written 

u^^ 2! 



u^ — 0^ oc 



a 



■y8 a-^e^ 



/33 a^ 



a^ 



•5 .y>3 



^^^' 



a 



•/9 3! 



r • • • 5 



or 



aV^* — /3V^' 



a 



/8 



Thus, we have 

Theorem 6. — 1% Theorem 4 the formula (31) wa;2/ ^^ replaced iy 



K(^) 



^2gaa: 



a^^"^ 






i^A^re a and /S are the roots of the quadratic x^ + KoX-i-Ki = 0. This form of 
K(x) is preferable for purposes of computation when kq and ki are of pre- 
ponderant importance as compared with k2, ks, .... 

We may in this way replace the terms of (31) to any extent by exponen- 
tials, and it is evident that when the terms thus replaced are those which 
involve kq, ki, ... tcn-i without involving k^, /c»+i, ..., the exponential terms 
which are obtained are precisely those which would be obtained as the value 
of K(a?) in Theorem 3, if k{x) were the polynomial 

/^.2 X^~^ 
fC0'^fCi(C-\'K2 — :+...4-/^»-i 7 rr-;. 

21 {n — l)\ 
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